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ABSTRACT 

We present Cosmic Origins Spectrograph (COS) measurements of metal abundances in eight 
0.083 < Zabs < 0.321 damped Lyman-a (DLA) and sub-damped Lya absorption systems serendipi- 
tously discovered in the COS-Halos survey. We find that these systems show a large range in metal- 
licities, with — 1.10< [Z/H] <0.31, similar to the spread found at higher redshifts. These low-rcdshift 
systems on average have subsolar metallicities, but do show a rise in metallicity over cosmic time 
when compared to higher-redshift systems. We find the average sub-DLA metallicity is higher than 
the average DLA metallicity at all redshifts. Nitrogen is underabundant with respect to a-group 
elements in all but perhaps one of the absorbers. In some cases, [N/a] is significantly below the lowest 
nitrogen measurements in nearby galaxies. Systems for which depletion patterns can be studied show 
little, if any, depletion, which is characteristic of Milky Way halo-type gas. We also identify affiliated 
galaxies for 3 of the sub-DLAs using spectra obtained from Keck/LRIS. None of these sub-DLAs arise 
in the stellar disks of luminous galaxies; instead, these absorbers may exist in galaxy halos at impact 
parameters ranging from 38 to 92 kpc. Multiple galaxies are present near two of the sub-DLAs, and 
galaxy interactions may play a role in the dispersal of the gas. Many of these low-redshift absorbers 
exhibit simple kinematics, but one sub-DLA has a complicated mix of at least 13 components spread 
over 150 km s~^. We find three galaxies near this sub-DLA, which also suggests that galaxy interac- 
tions roil the gas. This study reinforces the view that DLAs have a variety of origins, and low-redshift 
studies are crucial for understanding absorber-galaxy connections. 
Subject headings: galaxies: abundances — galaxies: ISM — quasars : absorption lines. 



1. INTRODUCTION 

Quasi-stellar object (Quasar or QSO) absorption lines 
have long been used to probe the intergalactic medium 
and also the interstellar medium (ISM) of galaxies at all 
redshifts. The damped Lyman-a systems (DLAs) (log 
A'hi > 20.3) are thought to contain mo st of the neutral 
gas in the Universe at high- r edshifts (IProchaska et al.l 
120091: INoterdaeme et al]|2009D . with sub-DLAs (19.0 < 
log A'hi < 20.3) also believed to hav e significant fraction 
at high-redshift (iPeroux et al.ll2005p. At lower redshift , 
a survey for sub-DLAs at z 2 bv lO'Meara et al.l ()2007[ ) 
found that they contain ^15% of the Universe's H atoms. 
The large H I column densities of both DLAs and sub- 
DLAs are typical of what is seen in the disks and nearby 
environments of galaxies. As such, these systems are be- 
lieved to be strongly related to galaxies at their observed 
redshift. At z ~ 0, D LAs have been associated to galac- 
tic ISM H waan et al.. ,2005i ). These systems provide a 
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unique means to probe galaxies over > 90% of the cos- 
mic history, without the biases introduced in magnitude- 
limited samples of galaxies studied in emission. 

Since these systems have such large column densities, 
the majority of the H is expected to remain neutral due 
to self-shielding. These high column densities also make 
absorption lines from metal atoms more detectable. Any 
metals inside one of these clouds will be dominantly in 
states with ionization potentials greater than that of H^ , 
since photoionization effects will be negligible. For sub- 
DLAs, the lower column densities make these photoion- 
ization effects more significant and corrections are typi- 
cally applied through use of ionization models. 

DLAs arc of particular importance to chemical evo- 
lution studies. Measurements of element abundances 
in such systems provide insights into the production of 
metals through cosmic time. This task is complicated 
by the large scatter in metallicities seen at any given 
redshift, due to the random sampling of galaxies with 
different masses and different environments. However, 
this is actually beneficial, since the sample will natu- 
rally probe a more representative range of all metallicities 
and provide an unbiased description of chemical enrich- 
ment. Despite the large scatter, using the large sample 
size of D LAs measured at high-redshift, IProchaska et al.l 
(|2003bf ) found evidence for an increase in elemental abun- 
dances with increasing cosmic time. This picture is con- 
sistent with our expectation that the Universe becomes 
more metal-rich as it evolves. Interestingly, many DLAs 
and sub-DLAs appear to be underabundan t in N with re- 
spect to the other metals in these systems (jJenkins et al.l 
[2005.: .Tripp et al...2005; .Pettini et al...2008l) . even at low- 
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rcdshifts. Observations of N in DLAs and sub-DLAs 
arc important for the study of N production and galac- 
tic chemical evo lution in gener al (iCenturion et al.ll2003l : 
iHenrv fc Procha ska 2007: Petti ni et al.ll2008| ). 

Unfortunately, relatively few measurements of 
low-redshift DLAs a. n d sub-DLAs hav e been made 
(jKulkarni et all 120051: iRao et all 120051) because the 
atmospheric cutoff of light below '^SOOOA makes ground 
based observing of the Lyman-a line impossible for red- 
shifts z < 1.65. With the Cosmic Origins Spectrograph 
(COS) aboard the Hubble Space Telescope (HST), 
we are able to study these low-redshift systems with 
unprecedented efficiency. This instrument provides us 
with a great opportunity to expand this small sample. 
These systems are crucial for a clear understanding of 
chemical evolution since z < 0.5 spans ^40% of the age 
of the Universe. 

Although much can be learned from absorption line 
data, the difficulty in the interpretation of these systems 
stems from the lack of information about the host galax- 
ies of DLAs and sub-DLAs at all rcdshifts. The number 
of DLAs with spectroscopically confirmed hosts is very 
small and this is in part due to the much lower number of 
low-redshift DLAs compared to high-redshift. At high- 
redshift, the cosmological dimming of the surface bright- 
ness of galaxies a (1 -|- z)~'^) makes galaxy detections 
more difficult. Low-redshift DLAs present a unique op- 
portunity to study the absorber-galaxy connection, since 
the host will be more easily detected. 

In this paper, we report on 3 DLAs and 5 sub-DLAs 
over a redshift range of 0.0830 < Zabs < 0.3211. In §2 we 
give background information on the discovery of these 
systems and on the observations. In §3 we describe our 
techniques for determining column densities from absorp- 
tion lines. In §4.1 we report the total column densities 
and abundances of each system, correction for photoion- 
izing effects in the sub-DLA cases. In §4.2 we report the 
depletion patterns of systems with a sufficient number 
of detected lines for analysis. In §4.3 the metallicities 
of our absorbers will be compared to those in the lit- 
erature and also with chemical evolution models. In §5 
we discuss nitrogen measurements for each system and 
possible implications on nitrogen nucleosynthesis. In §6 
wc identify affiliated galaxies for 3 of the sub-DLAs us- 
ing observations from a redshift survey with Kcck/LRIS. 
Conclusions arc summarized in §7. Throughout this pa- 
per, we adopt a cosmological model with fi^ — 0.30, 
Qa = 0.70, and Ha = 70 km s^^ Mpc-\ 

2. PROGRAM SUMMARY AND OBSERVATIONS 

The DLAs and sub-DLAs in this sample were dis- 
covered serendipitously as part of programs G011598 
and G012248 (PI: Tumlinson), programs focused on 
studying multiphase baryons in halos oi L > L* galax- 
ies and dwarf galaxy halos, respect ively. Early result s 
from these progr a ms can be fo und inlThom "eFal|(l20Tl . 
iTumlinson et al.l (|2011aD and iTumlinson et all (1201111 ). 
The names and basic properties of these absorbers are 
presented in Table [T] The DLA with the highest col- 
umn density in our sample, SDSS J1009-I-G713, was pre- 
sented in detail in iMeiring et al.l ()2011[ ) along with ba- 
sic information of the other systems. The details of our 
COS observation s are s vmimarized in Table [H Follow- 
up Keck/HIRES (jVogt ct al...l994 ) observations were ac- 



Table 1 

Properties of the DLAs and sub-DLAs 



QSO 




'^abs 


log A^Hi 








(cm~^) 



SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 



J092554. 
J092837. 
J100102. 
J100902. 
J143511. 
J155304. 
J161649. 
J161916. 



70-1-400414.1 
98-1-602521.0 
55+594414.3 
06-1-071343.8 

53- 1-360437.2 
92-1-354828.6 
42-1-415416.3 

54- 1-334238.4 



0.471 
0.295 
0.746 
0.456 
0.429 
0.722 
0.440 
0.471 



0.2477 
0.1538 
0.3035 
0.1140 
0.2026 
0.0830 
0.3211 
0.0963 



19.55±0.15 
19.35±0.15 
19.32±0.10 
20.68±0.10 
19.80±0.10 
19.55±0.15 
20.60±0.20 
20.55±0.10 



Note. Exeluding JlOOl-l-5944, measurements 
et al. (2011). Second and third columns 
of QSO emission and the redshift of the 
respectively. 



taken from Meiring 
denote the redshift 
absorption system, 



Table 2 

Summary of COS and Keck observations 



QSO 


Obs. Date 


Grating 


Central A 


Exp. Time 








(A 


) 


(s) 


J0925-I-4004 


2010-06-15 


G130M 


1318, 


1327 


2176, 1589 




2010-06-15 


G160M 


1600, 


1623 


1248, 3055 


J0928+6025 


2010-03-06 


G130M 


1291, 


1309 


1163, 1148 




2010-03-06 


G160M 


1600, 


1623 


1547, 1505 


JlOOl-l-5944 


2011-03-18 


G130M 


1291, 


1309 


1618, 1682 




2011-03-19 


G160M 


1600, 


1623 


2710, 2491 


J1009-I-0713 


2010-03-29 


G130M 


1291, 


1309 


2191, 1497 




2010-03-30 


G160M 


1600, 


1623 


2007, 2002 


J1435+3604 


2010-05-22 


G130M 


1291, 


1327 


1788, 1678 




2010-05-22 


G160M 


1577, 


1623 


2212, 2212 


J1553-I-3548 


2010-10-01 


G130M 


1291, 


1309 


1600, 1623 




2010-10-02 


G160M 


1600, 


1623 


1436, 1401 


J1616-I-4154 


2010-07-19 


G130M 


1291 


3644 




2010-07-19 


G160M 


1600, 


1623 


2437, 2167 


J1619-I-3342 


2010-08-19 


G130M 


1291, 


1309 


2358, 2989 




2010-08-19 


G160M 


1600, 


1623 


4394, 4395 


J0928+6025 


2010-03-26 


HIRES 






3000 


J1009+0713 


2010-03-26 


HIRES 






3600 


J1619-I-3342 


2010-09-02 


HIRES 






3600 



quired for three of our absorption systems to see more 
detailed component structure, since COS has lower reso- 
lution (R~20,000) compared to HIRES (R-44,000), and 
those are also summarized in table [5] A more detailed 
description of the COS observations for this project 
along with its data reduction methods were outlined in 
IMeiring et all (poTll ). 

After the data reduction, the spectra were binned by 
3 pixels, as the raw COS data are oversampled with a 
~6 pixel wide resolution element. All subsequent mea- 
surements and analysis were performed on the binned 
spectra. The binned spectra have a resolution of ^15 
km s~^ per resolution element. The signal-to-noise ra- 
tio (S/N) of these COS spectra range from 7 to 15 per 
resolution element. 

Accurate column densities for the H I A 1215.67 ab- 
sorption line in DLAs and sub-DLAs are measured by 
fitting the damping wings in the profile. Fits to the Ly- 
a line s all of our systems were presented in lMeiring et al.l 
(|2011f) . with the exception of J1001-(-5944 which is pre- 
sented in Figure[TJ The Ly-a lines for all systems, except 
JlOOl-l-5944, were fitted using a single component with 
the continuum and profile fit simultaneously. The red- 
shifts of the systems were determined from the metal 
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J1001 + ;5944 log N(H I) = 1 9,32 ±0, 10 z„|„=0,3035 
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Figure 1. Fit to the Ly-a line for JlOOl+5944. The sohd green 
line represents the adopted continuum and the red hnes are the 
best-fit model with column densities la above and below the best- 
fit value. The Icr flux uncertainty is shown in cyan below the data. 

lines, which with narrow features provide for a more ac- 
curate redshift determination. For JlOOl+5944, we used 
the four components found in O I Hnes to fit the Ly-a 
hne. In total, we find 3 DLA systems and 5 sub-DLA 
systems in these COS data. Absorption rcdshifts and 
neutral hydrogen column densities for the systems were 
given in Table [1] Errors on A^hi measurements given in 
Table [1] have been estimated by eye and are dominated 
by the continuum placement uncertainty. 

3. ANALYSIS 

Here we present a detailed description of our meth- 
ods for determining the column densities from observed 
absorption features. The column densities for unsatu- 
rated lines were estimated from a multicomponent Voigt 
profile fit using the fitting code of iFitzpatrick fc Spitzei] 
(|1994f ). This code minimizes the value between the 
data and the synthetic Voigt profiles convolved with the 
COS line spread functions (LSPfl at the nearest tabu- 
lated wavelength (jChavamian et al1l2009() . The column 
density errors for our profile fits are formal Icr uncer- 
tainties from the minimization and do not include 
contributions from the uncertainties in the continuum 
placement. The fitting procedure for each system was 
to select the species with the most extensive number of 
observed lines, preferably over a large range of oscilla- 
tor strengths, /, to reduce the possibility of saturation 
effects. Then we use the velocity ccntroids and Doppler 
parameters, b, from the best fit of this species for the 
fits to other lines and only allow the column densities 
in each component to vary (see appendix for more de- 
tails). We fit multiple transitions simultaneously when- 
ever possible. Lines showing evidence of being strongly 
blending with another line were not used for the profile fit 
of that transition. The COS spectrum was also allowed 
to shift in velocity space to account for calibration issues 
in the instrument, which can cause offsets between var i- 
ous transitions of order ±10km (jSavage et al.ll201ll ). 
If Keck data was available, then we determined the com- 
ponent structure by fitting the Mg I A 2852.96 line. The 
Keck/HIRES spectra were convolved with a Gaussian 
LSF with FWHM= 2.5 pixels. 

It is worth noting that the limited resolution of 

Tabulated LSFs for COS are available at 
http : / /www. stsci . edu/hst/cos/perf ormance/spectral.resolutlon/ 



COS will prevent it from detecting components that 
Keck/HIRES can resolve. This will not dramatically 
affect column density measurements but does hinder 
our understanding of the kinematics of these systems. 
Component fits using only the COS spectrum should be 
treated with caution for analysis of kinematic structure 
as they likely do not refiect the true number of compo- 
nents. 

Our profile fitting approach is chosen over apparent op- 
tical depth (AOD) methods whenever possible because 
the AOD is more prone to systematic errors from mild 
saturation and noise; both of which are concerns in these 
COS spectra due to the relatively low S/N and broad LSF 
of the COS spectrograph. By fitting multiple transitions, 
over a range of / values, we have a better chance to over- 
come these effects. We assume that species with similar 
ionization potentials can be fit using identical b values be- 
cause the resolution of COS is not high enough for small 
changes in b, caused by the different atomic masses, to 
be significant. One also expects different species with 
similar ionization potentials to exist in similar regions 
of a cloud and thus have similar velocity components. 
Intermediate-ions, such as Fe^^ and Si^^, appear to 
have velocity structure that is indistinguishable from the 
low-ion profiles in most systems and so it has been argued 
that the two regions are physically connected or are par- 
tially mixed (jPessauges-Zavadskv et aT]|2003D . For this 
reason we also choose to fit the intermediate-ions using 
the same parameters from the fit to the lower ions. Rest 
fram e wavelengths and / values for all lines were taken 
from lMortonI |2003[). 

Lines were determined to be saturated in two ways: (1) 
We compared the total column density found from pro- 
file fitting with the apparent column density found from 
direct integration of the line using the AOD method of 
iSavage fc Sembachl (|1991[ ). Cases for which apparent col- 
umn density measurements and profile fitting measure- 
ments were inconsistent within their uncertainties were 
determined to be saturated. (2) For species with several 
lines detected over a range of / values, lines with larger 
/A values will give lower values of apparent column den- 
sity than lines with smaller /A values if saturation occurs. 
When a line is determined to be saturated we report a 
lower limit based on the apparent column density found 
for the line of that species with the lowest /A value (i.e. 
the least saturated line). 

Since the component structure in the more highly ion- 
ized gas traced by the CIV AA 1548.20, 1550.77, Si IV 
AA 1393.75, 1402.77 and O VI AA 1031.92, 1037.61 lines 
are often different than that of the low ions, we have 
measured the column density of these species by direct 
integration of the line using the apparent column density. 

We present the adopted velocity components and b val- 
ues for each system as well as the column densities of 
detected species at each of these component in the ap- 
pendix. The corresponding fits using these parameters 
are shown in Figure [2lfT0l 

4. CHEMICAL ABUNDANCES 

4.1. Total Column Densities and Abundances 

Total column densities measurements were determined 
by summing the individual components from Table [TT] 
or from apparent column density measurements and are 
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Table 3 

Total logarithmic column densities, log A^x (cm^^), for each absorber 



X 


J0925+4004 


J0928+6025° 


JlOOl+5944 


J1009+0713'' 


J1435+3604 


J1553+3548 


J1616+4154 


J1619+3342 


HI 


19.55±0.15 


19.35±0.15 


19.32±0.10 


20.68±0.10 


19.80±0.10 


19.55±0.15 


20.60±0.20 


20.55±0.10 


CII 


>15.18 


>14.91 


>15.06 


>15.90 


>14.52 


>14.35 


>15.15 


>14.38 


cm 


>14.17 




>14.52 




>14.33 








CIV= 




>14.1 




13.86±0.20 




13.99±0.08 




13.28±0.08 


CII* 


<14.18 


<13.82 


<13.73 


<13.70 


<13.68 


<13.46 


13.95±0.08 


<13.12 


NI 


14.75±0.04 


14.10±0.09 


13.65±0.17 


15.11±0.10 


14.60±0.14 


<13.74 


>15.35 


14.64±0.12 


Nil 


>14.95 


>14.85 


>14.95 






14.16±0.07 


>15.17 


>13.92 


OI 


15.95±0.09 


>15.08 


15.64±0.02 


>16.00 


>15.58 


>14.56 


>15.45 


>14.80 


OVP 






14.34±0.05 












Mgl 




12.70±0.06 




>12.85 








12.40±0.14 


Mgll 




>13.99 




>15.35,<16.0 










Sill 


14.62±0.06 


>14.39 


14.73±0.03 


>15.00 


>14.11 


14.22±0.05 


>15.08 


>13.93 


Silll 


>13.74 


>13.77 


>14.00 




>13.41 


>13.3 


>13.93 


>13.24 


SilV^ 


13.54±0.11 


13.86±0.12 






13.20±0.13 


13.59±0.08 




13.06±0.05 


PII 


<13.50 


<14.17 


12.81±0.07 


<13.80 


<12.79 




13.46±0.10 


13.16±0.19 


SII 


<14.72 


<14.65 


<14.53 


15.25±0.12 


14.60±0.12 


<14.24 


15.37±0.11 


15.08±0.09 


Call 




12.81±0.06 




11.96±0.08 








12.42±0.02 


Till 




<11.94 




12.70±0.03 








11.90±0.04 


Fell 


14.22±0.09 


14.90±0.08 


14.30±0.04 


15.29±0.17 


14.20±0.08 


14.01±0.07 


15.02±0.05 


14.38±0.15 


Felll 


<14.23 


14.59±0.11 


14.14±0.09 




<13.69 




14.51±0.04 


13.95±0.21 


Nill 


<14.41 


<13.67 


<14.23 


13.93±0.18 


<13.99 


<13.92 




<13.53 



soft UV re gion for di s tance s as far as 100 kpc from the 
galaxy (see lFox et all (|2005[ ) for further details). CLOUDY 
models were also made using just the estimated extra- 
galactic background flux from iHaardt fc Madaul (|19960 
to determine if there would be significant differences. 
Even with the different shapes in the flux curves for 
these sources, the two models give almost identical re- 
sults. Therefore, we feel it is reasonable to continue with 
our original approach using the galactic radiation. 

Oxygen is the best indicator of metallicity in sub-DLA 
systems because it traces the H° region due to charge 
exchange reactions between O and H. For this reason, its 
abundance was used as the initial estimate for metallicity 
in the CLOUDY models. If accurate O measurements were 
unavailable, we used S or Si. For the case of J0928-I-6025 
we were forced to use Fe. Since Fe is prone to deple- 
tion, we will apply additional corrections to metallicity 
after modeling. Using these mctallicities, CLOUDY mod- 
els of the total column densities of elements as a func- 
tion of ionization parameter, [/, were made. The allowed 
range of U is constrained using the measured values of 
the column densities for observed elements along with 
the observed ratio of Fe++/Fe+ and Si++/Si+. If the 
intermediate-ions are physically connected with the low- 
ions, as was discussed earlier, then these ratios should 
constrain the value of U and also be less affected by dust 
depletion. An extensive study on ioni zation corrections 
in a m ultiphase ISM using CLOUDY by iMilutinovic et ahl 
(|20T0l) found that using the A1'*"'"/A1+ ratio recovers the 
correct log U value to within 0.1 dex for a two-phase 
medium. 

For the majority of cases, the species N°, Si+, Ti+, and 
Fe"*" appear to be underabundant relative to predicted 
CLOUDY abundances. The underabundance of Si, Ti, and 
Fe is likely due to their depletion onto dust grains. The 
underabundance of N will be discussed in the next sec- 
tion and is likely the result of its nucleosynthetic origins. 



The COS lines were determined using the apparent column density. 
Measurements taken from Meiring et al. (2011). 
Higher ions measured using apparent column density. 

listed in Table [31 Upper limits to undetected lines were 
derived from a Str upper limit on the equivalent width, 
assuming the line is on the linear part of the curve of 
growth. 

The absolute abundances were measured with respect 
to solar using the following formula; 

[X/H] = log[N(X)/N(H)]DLA - [X/H]q + ex (1) 

where ex is an ionization correction factor applied to the 
sub-DLAs. Sola r system reference ab undances ([X/H]©) 
were taken from lAsplund et all ()2009D adopting the pho- 
tospheric values for C, N and O and meteoric values for 
all other elements. In these systems the dominant species 
is the neutral state for O and N and the first ions for C, 
Mg, Si, P, S, Ti, Fe, and Ni. Abundance measurements 
are not made for Ca, since the dominant ion in this case is 
Ca"'"^, which doesn't have lines in our observing window. 
For saturated and unobserved lines we present lower and 
upper limits to abundances, respectively, using the AOD 
method. In Table U] we show the resulting abundances 
for each absorber without applying ex. 

Since the column densities of sub-DLAs are not high 
enough to keep the H gas mainly neutral through self- 
shielding, it is necessary to correct for effects of pho- 
toionization through models. It has been suggested that 
corrections sh o uld be small fo r A'hi > 19.5 at low-redshift 
(jViegasI [19951 : [Jenkinj I2009D . However, for complete- 
ness, we estimated corrections for all of the sub-DLAs. 
Ionization corrections for sub-DLA systems are deter- 
mined using the CLO UDY ionization code (version 08.00, 
iFerland et al.l (jl998D '). In this model, the gas is consid- 
ered to be a uniform slab that is illuminated from one side 
by photoionizing radiation. We adopt the ionizing flux 
in the disk o f the Milky Way based on the calculations 
of iFox et all (l2()05) . If sub-DLAs are associated with 
nearby galaxies, the fiux from starlight in these galaxies 
should dominate over the cxtragalactic background in the 
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Table 4 

Abundances, [X/H], for each absorber with respect to solar 



X 


J0925+4004'' 


J0928+6025° 


JlOOl+5944'' 


J1009+0713'' 


J1435+3604'' 


J1553+3548" 


J1616+4154 


J1619+3342 


c 


>-0.80 


>-0.87 


>-0.69 


>-1.2 


>-1.72 


>-1.63 


>-1.83 


>-2.60 


N 


-0.63±0.16 


-1.08±0.17 


-1.50±0.20 


-1.40±0.14 


-1.03±0.17 


<-1.64 


>-1.08 


-1.74±0.16 


o 


-0.29±0.17 


>-0.96 


-0.37±0.10 


>-1.40 


>-0.91 


>-1.68 


>-1.84 


>-2.44 


Mg 




>-0.89 




>-0.86,<-0.21 










Si 


-0.44±0.16 


>-0.47 


-O.lOiO.lO 


>-1.2 


>-1.20 


-0.84±0.16 


>-1.03 


>-2.13 


P 


<0.52 


<1.39 


0.06±0.12 


<-0.31 


<-0.44 




-0.57±0.22 


-0.81±0.21 


s 


<0.02 


<0.15 


<0.06 


-0.58±0.16 


-0.35±0.16 


<-0.46 


-0.38±0.23 


-0.62±0.13 


Ti 




<-0.32 




-0.89±0.10 








-1.56±0.11 


Fe 


-0.90±0.19 


0.10±0.17 


-0.47±0.11 


-0.84±0.20 


-0.83±0.15 


-0.99±0.17 


-1.03±0.21 


-1.62±0.18 


Ni 


<0.66 


<0.12 


<0.71 


-0.95±0.20 


<-0.01 


<0.17 




<-1.22 



" Without applying photoionization correction factors. 
^ Measurements taken from IMeiring et al.l II2011I ). 



Table 5 

Ionization corrections, exi for the sub-DLA absorbers 



X 


J0925+4004 


J0928-I-6025 


JlOOl+5944 


J1435-f-3604 


J1553-I-3548 




-3.4<log U<-2.S 


-3.9<log U<-3.5 


-3.6<log U<-3.3 


log U<-3.3 


-3.8<log U<-2.7 


c 


-0.40,-0.22 


-0.22,-0.12 


-0.37,-0.27 


-0.16, 0.03 


-0.48,-0.15 


N 


0.01, 0.03 


0.01, 0.02 


0.02, 0.03 


0.00, 0.07 


0.01, 0.04 


Mg 


-0.30,-0.17 


-0.19,-0.10 


-0.30,-0.21 


-0.12, 0.04 


-0.33,-0.11 


Si 


-0.45,-0.25 


-0.26,-0.14 


-0.39,-0.28 


-0.18, 0.00 


-0.38,-0.14 


P 


-0.37,-0.21 


-0.21,-0.12 


-0.37,-0.27 


-0.15,-0.01 


-0.48,-0.16 


s 


-0.25,-0.16 


-0.17,-0.11 


-0.28,-0.21 


-0.11, 0.00 


-0.30,-0.12 


Ti 


-0.11,-0.10 


-0.13,-0.08 


-0.18,-0.15 


-0.06,-0.01 


-0.08,-0.07 


Fe 


-0.16,-0.11 


-0.13,-0.08 


-0.19,-0.16 


-0.07, 0.00 


-0.14,-0.08 


Ni 


-0.28,-0.21 


-0.22,-0.13 


-0.35,-0.28 


-0.15, 0.00 


-0.32,-0.17 



Note. Values shown are the limits for ex over the adopted range of log U. 



Table 6 

Adopted Metallicity for each absorber 



Absorber 


^abs 


log Nui 


Z 


[Z/H] 






(cm~^) 







J0925-I-4004 


0.2477 


19.55±0.15 


O 


-0.29 ±0.17 


J0928+6025 


0.1538 


19.35±0.15 


Fe 


-f-0.31 ± 0.17"'' 


JlOOl+5944 


0.3035 


19.32±0.10 


O 


-0.37 ±0.10 


J1009+0713 


0.1140 


20.68±0.10 


s 


-0.58 ±0.16 


J1435+3604 


0.2026 


19.80±0.10 


s 


-0.41 ± 0.16" 


J1553+3548 


0.0830 


19.55±0.15 


Si 


-1.10 ±0.16" 


J1616+4154 


0.3211 


20.60±0.20 


s 


-0.38 ±0.23 


J1619+3342 


0.0963 


20.55±0.10 


s 


-0.62 ±0.13 



" After applying average ex for element used. 

^ Applie d a +0.3 dex correction, as discussed in IHafelski et al.l 
din prepi ). 



It is also apparent that the measured abundances of the 
higher ions, such as Si'^"'" and C'^"'", do not appear to agree 
well with the lower ions in the CLOUDY models. How- 
ever, as was mentioned earlier, the higher ions may be 
located in a different region than that of the lower ions 
and it would be unwise to use these states as a guide in 
determining the conditions of the lower ions. For sev- 
eral of our systems the kinematics of these high-ions do 
show significant difference from the low-ions, suggesting 
a multiphase gas component. The high-ions are still use- 
ful because the value of U in those regions will certainly 
be higher than the low-ion regions and can be used as a 
conservative upper limit to the possible ionization. 

For J0925+4004, we used the upper limit of Fe++ /Fc+ 
and the lower limit of Si+^/Si^ to constrain log U 
and this also agrees with the limit imposed from Si'^+. 



For J0928+6025, JlOOl+5944, J1435+3604 we use the 
Fe++/Fe+ ratio to constrain log U. For JlOOl+5944, 
this region agrees with the observed column densities of 
Fe++, Si+ and P+. For J1553+3548, we used the up- 
per limit suggested by the Si'^^ column density and the 
lower limit of Si+^/Si^ to constrain U. The range of 
values we adopt for log U can be seen in Table [5l along 
with the corresponding range in ex- We will adopt the 
average value of ex for each individual element when de- 
termining the metallicity for these systems. For exam- 
ple, in J1553+3548 we find -0.38 < esi < -0.14 and so 
the correction we will use for this system's metallicity, 
which is determined from Si in this case, is esi = —0.26. 
Since O has negligible ionization correction because of 
its relationship to H, the metallicity measurements for 
sub-DLAs J0925+4004 and JlOOl+5944 are independent 
from the degree of photoionization. 

The metallicities adopted for the DLAs J1616+4154 
and J1619+3342 were based on the abundance derived 
from the S II lines because this ion provides accurate 
column density measurements and it has relatively low 
depletion rate (however, see the caveats discussed by 
iJenkinsI ()2009f )V For both cases, this value can be com- 
pared to the abundance found for P, which also has low 
depletion levels. For J1616+4154 and J1619+3342, we 
found [S/P]=0.19±0.16 and [S/P]=0.19 ± 0.21, respec- 
tively. The metallicities found from either of these ele- 
ments appears consistent with the other within the un- 
certainties. However, it may suggest a enhancement. 
The metallicities we adopt for each of our absorbers is 
shown in Table [B] 
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Figure 2. Velocity plots for several lines in the COS spectra of J0925+4004. The dashed lines indicate the velocity components, the data 
is presented in black, the error array is given in blue, and the best fit is brown. The Fe II A 1144.94 hne required a large offset to fit (—40 
km s~^) suggesting that it might be blended with something. The following lines were not used for profile fits: N I AA 1200.22 and 1200.71, 
since they are significantly blended together; Si II A 989.87 appears blended, perhaps with N III A 989.80. Unrelated features or blends are 
indicated with magenta circles. Rest velocity corresponds to 2: = 0.2477. 
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Figure 3. Velocity plots for several lines in the HIRES spectra of 
J0928+6025. The dashed lines indicate the velocity components, 
the data is presented in black, the error array is given in blue, and 
the best fit is brown. Rest velocity corresponds to z = 0.1538. 

4.2. Depletion Patterns 

The depletion patterns for these systems were esti- 
mated using the model of depletion pa tterns in the ISM 
of the Milky Way from lJenkiiii (|2009f ). In this model a 
single value (F*) indicates the overall level of depletion 
in a sightline, with larger values of indicating higher 
levels of dust depletion. Accurate determination of F* 
requires a large number of elements which are spread 
over a broad range of condensation temperatures. Only 
a fraction of these lines are available to us because of the 
low S/N of our data and limited wavelength coverage. 
Despite this, it is still possible to see general trends the 
depletion patterns of these systems. 

We are able to estimate a depletion factor of 
F* for J0925+4004, JlOOl+5944, J1435+3548 and 
J1619-I-3342, as these systems have a sufficient num- 
ber of lines for basic analysis. For J0925-I-4004 we find 
F* = —0.13 ± 0.19, implying an intrinsic metallicity of 
[Z/II]=: —0.41 ± 0.15, which is consistent with the metal- 
licity determined from the O I lines, [0/II]= -0.29±0.17. 
For J1001-h5944 we find F* = -0.69 ± 0.15, implying 
an intrinsic metallicity of [Z/II]=: —0.79 ± 0.13, which is 
below the value determined from the O I line, [0/H] = 
-0.37±0.10. The depletion pattern for JlOOH-5944 does 
not seem to fit the linear trend of F* very well, since there 



is a small number of diagnostic lines, this is likely the 
source of discrepancy. The deficiency of N in this system 
may also contribute to the problem. For J1435-I-3548 
we find F* = —0.78 ± 0.22, implying an intrinsic metal- 
licity of [Z/H]= -0.97 ± 0.21, which is below the value 
determined from the S II line, [S/H]= -0.41 ± 0.16. For 
J1619+3342 we find F* = -0.55 ± 0.11, implying an in- 
trinsic metallicity of [Z/II]= —1.53 ± 0.17, which is well 
below the metallicity determined from the S II and P II 
lines, [S/II]= -0.62 ± 0.13. For the latter two systems, 
one possible reason for the discrepancy is that both P 
and S show some what anomalou s behavior in the way 
they deplete (i.e. iJenkinsI (120090 ). Additionally, some 
of the S II could arise from regions where the hydrogen 
is ionized, since the ionization potential of S II is 23.4 
eV, which could create a false impression that the abun- 
dance of S relative to H is higher than reality. All of 
our systems show a trend towards negative values of F.^ 
which indicates little, if any, depletion. This is typically 
seen along sightlines with extr emely low ave rage densi- 
ties in the halo of our Galaxy ()Jenkinj|2009() . implying 
that each of these sightlines pass through the outer edge 
of a galaxy. 

4.3. Metallicity Evolution of DLAs and Sub-DLAs 

The absorption systems in this study show a range of 
metallicities from above solar for J0928-t-0625, to one- 
tenth solar for J1553+3548. We now compare the ob- 
served metallicities of our low-redshift DLAs and sub- 
DLAs to those in the literature. The nearly undepleted 
clement Zn is a commonly adopted metallicity tracer 
in DLA and sub-DLA systems and is used for some of 
the measurements in the sample we constructed (How- 
ever, see see the caveats noted by iHowk fc SembachI 
(|1999f )). Most of the other measurements are based 
on abundances of a-elements (O, Si, S). The use of a- 
elements with Zn measurements does create an inhomo- 
geneous sample because of the different nucleosynthctic 
origin s between these elements, however iKulkarni et al.l 
(|2007f ) found that their inclusion in the samples made 
essentially no change in their results. The findings by 
IProchaska fc Wolfd (l200l ) that [Si/Zn]= 0.03 ±0.05 also 
suggests that combining these populations seems reason- 
able. The remaining sample is made up of metallicities 
determined from Fe using a -1-0.3 dex correction for de- 
pletion. 

The comparison sample we will utilize includes 194 
DLAs with 0.0096 < z < 4.743 and 53 sub-DLAs with 
0.0063 < z < 3.1710. The majo rity of DLA so u rces com e 
from the list constructed by IRafelski et al.l (jin prepO . 
This list is designed to be unbiased at z > 1.5. The 
major ity of sub-DLA sources co mes from IKulkarni et al.l 
(|2007f) and lMeiring et al.l (|2009D . The DLA sample was 
divided into 10 redshift bins with equal numbers of 
sources in each bin. The unweighted and weighted mean 
metallicity was determined for each bin and is shown in 
Table [71 Weighted metallicities are defined as: 



^w^[Z/H], 



(2) 



Where the weight, Wi, is the reciprocal square of the 
corresponding uncertainty (wi = 1/af). This was re- 
peated for the sub-DLAs over 4 redshift bins and these 
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Figure 4. Velocity plots for several lines in the COS spectra of J0928+6025. The dashed lines indicate the velocity components, the data 
is presented in black, the error array is given in blue, and the best fit is brown. Unrelated features or blends are indicated with magenta 
circles. O I A 1039.23 is blended with N I A 1199.55 in the MW. Si IV A 1393.76 is blended with Fe II A 1608.45 line in the Milky Way 
(MW). Rest velocity corresponds to 2: = 0.1538. Fits were not made for this system because of the low S/N. 



are shown in Table [H The unweighted binned values 
can be seen alongside chemical evolution models from 
Dave &: Oppenheimeii (j2007D (magenta dashed line) and 



Pei et all (|l999l ) (blue dash-dot line) in Figure HH The 



look-back times are estimated assuming a cosmological 
model with = 0.30, fl/^ = 0.70, and Hq ^ 70 
km Mpc^^. For each bin, the data points arc plot- 
ted at the median value of the bin and the horizontal 
error bars denote the range in look-back time spanned 
by that bin. The vertical error bars denote the uncer- 
tainty in the mean metallicity and accounts for scatter 
and also measuremen t uncertainties. The model from 
iDave fc Oppenheimeii (|2007[ ) is a hydrodynami c simu- 
lation whereas the mod el from iPei et al.l (|1999[ ) model 
is semi-analytical. The IPei et al.l (|1999f ) model agrees 
with the DLA systems at Zabs > 1-8, but overpredicts 
the mean metallicity of low-redshi ft systems. The model 
from IDave fc Qppenheimerl ()2007t) only agrees with low- 
redshift DLA observations, and seems to overpredict the 
mean mctallicities in the highest rcdshift systems. The 
binned sub-DLA data agree with both models, however 
there is a deficiency in data points at high-redshifts. 

As expected, the low-redshift systems do appear to 
have higher average mctallicities than higher redshift 



systems. We also see a relatively large scatter in 
metallicity present among our sources as we would 
expect from random probing of galaxies without re- 
gard to mass. As was the case in previo us studies 
(jKulkarni et"all l2007t iProchaska et al.ll2003bl ). we con- 
tinue to see an underabundance in mean metallicity at 
low-rcdshift relative to solar. This trend shows agree- 
ment with recent cosmological hydrodynamic simulations 
(jpavc & Oppcnhcimcr 2007; Cen 2 01 Ol ). It is also ap- 
parent that sub-DLAs, on average, show slightly higher 
mean mctallic ities than DLA s , but not to extent found 
previously bv iKulkarni et al.l ()2007[) . However, the rea- 
sons for this trend are still unclear. 

5. NITROGEN NUCLEOSYNTHESIS 

The nucleosynthetic origins of nitrogen is a topic that 
has proved to be quite complicated. The production of 
N occurs from primary and secondary production pro- 
cesses in intermediate mass stars. Primary production 
occurs from newly synthesized C via the CNO cycle 
and secondary production uses C and O f rom p revious 
star formation in the CNO cycle (jMarigol [20011) . Mea- 
surements of high-redshift DLAs, where N I lines are 
redshifted into wavelengths where ground-based observ- 
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Figure 5. Velocity plots for several lines in the COS spectra of JlOOl+5944. The dashed lines indicate the velocity components, the data 
is presented in black, the error array is given in blue, and the best fit is brown. Unrelated features or blends are indicated with magenta 
circles. The following lines were not used for profile fits: O I A 971.74 is blended with the hyy line of the same system; O I A 988.77 appears 
to be blended with and unidentified line based on its strength. Rest velocity corresponds to z = 0.3035. 
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Figure 6. Velocity plots for several lines in the COS spectra of J1435+3604. The dashed lines indicate the velocity components, the data 
is presented in black, the error array is given in blue, and the best fit is brown. Unrelated features or blends are indicated with magenta 
circles. Si II A 1020.70 is blended with the Ly7 line of an absorber at 2:^65 = 0.2624. Si IV A 1402.77 is blended with the Lya line of an 
absorber at z^bs = 0.3874. Rest velocity corresponds to z = 0.2027. 
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Figure 7. Velocity plots for several lines in the COS spectra of J1553+3548. The dashed lines indicate the velocity components, the data 
is presented in black, the error array is given in blue, and the best fit is brown. Unrelated features or blends are indicated with magenta 
circles. The Si IV A 1393.76 appears to be blended with an unidentified line. Rest velocity corresponds to z = 0.0830. 



Table 7 Table 9 

Binned DLA Mean Metallicities [N/a] and [a/H] for each absorber 



Redshift Bin 


# of systems 


Unweighted 


Weighted 


ID 


a 


[N/a] 


[a/H] 


0.0096-0.7405 


19 


-0.73±0.46 


-0.68±0.46 


J0925-1-4004 


O 


-0.32±0.10" 


-0.29±0.17 


0.7731-1.8644 


19 


-0.99±0.33 


-1.02±0.33 


J0928-I-6025 


O 


<-0.10" 


>-0.67 


1.8918-2.1408 


19 


-1.07±0.50 


-1.13±0.50 


JlOOl-l-5944 


o 


-1.10±0.17" 


-0.37±0.10 


2.1410-2.3700 


19 


-1.34±0.47 


-1.20±0.47 


J1009-I-0713 


s 


-0.82±0.16 


-0.58±0.16 


2.3745-2.5841 


19 


-1.45±0.58 


-1.18±0.58 


J1435-I-3604 


s 


-0.58±0.18'' 


-0.41±0.16" 


2.5950-2.7835 


19 


-1.21±0.49 


-1.23±0.49 


J1553-I-3548 


Si 


<-0.51" 


-1.10^10.16" 


2.7958-3.0621 


20 


-1.44±0.50 


-1.26±0.50 


J1616-I-4154 


s 


>-0.70 


-0.38±0.23 


3.1040-3.3655 


20 


-1.59±0.56 


-1.57±0.56 


J1619-I-3342 


s 


-1.12±0.15 


-0.62±0.13 


3.3869-3.7609 


20 


-1.69±0.52 


-1.75±0.52 


" Applied average ex for elements used. 



3.7612-4.7430 20 -1.78±0.48 -1.78±0.48 



Table 8 

Binned sub-DLA Mean Metallicities 



Redshift Bin 


# of systems 


Unweighted 


Weighted 


0.0063-0.7821 


13 


-0.20±0.60 


-0.36±0.60 


0.8426-1.1157 


13 


-0.41±0.59 


-0.51±0.59 


1.1414-1.4094 


13 


-0.55±0.37 


-0.64±0.37 


1.4244-3.1710 


14 


-0.80±0.73 


-0.43±0.73 



ing is possib le, find N to be underabundant relative to 
g-elements (ICenturion et aLl 119981 : iHenrv fc Prochas^ 
I2007t iPettini et all 120081) . The study of N relative to 
a-process elements is useful because it can give insights 
into the star formation history of a gas system, since 
[N/a] will increase as a system ages and intermediate 



mass stars have time to inject more N into the system. 

The a-process elements that we will be utilizing are O, 
Si, and S. The specific a-element used was determined 
for each system in the following way. The primary choice 
was O, because of its resonant exchange with H. The next 
choice for a was S, since it shows only mild depletion onto 
dust grains. If neither O or S are accurately determined, 
then the abundance for Si is adopted. However, since Si 
shows a strong tendency to deplete onto dust grains this 
estimate is less certain. The measured values for [N/a] 
arc listed in Table [HI It appears to be the case for the 
majority of our absorption systems that N is underabun- 
dant relative to a-process elements. 

Abundance measurements of N have been measured 
previously using optical spectra of H II regions of nearby 
low- luminosity galaxies, whe re N is mainl y prod uced 
from primary production, by IVan Zee et al.l ()2006l) and 
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Figure 8. Velocity plots for several lines in the COS spectra of J1616+4154. The dashed lines indicate the velocity components, the data 
is presented in black, the error array is given in blue, and the best fit is brown. Unrelated features or blends are indicated with magenta 
circles. The N I lines show strong blending with each other. SUA 1253.81 is blended with C I A 1656.93 in the MW. SUA 1259.52 is 
blended with Si II A 1260.42. Rest velocity corresponds to z = 0.3213. 



show a flat plateau in [N/0] ratio with respect to metal- 
licity. This ratio was also measured in higher metal- 
licity H II regions in nearby spiral g alaxies, where 
more secondary production occurs, by IVan Zee et alj 
(jl998| ) and these show a steady increase in [N/0] for 
larger metallicities. In Figure [121 we show [N/a] 
vs. [a/H] for our sample along with these H II re- 
gions, as we ll as t he high-redshift DLA sample from 
iPettini et all (l2008f). th e low-redshift sub-DLA system 



from Tripp et alJ (120051) . and a lo w-redshift DLA from 
iBowen et all (|2005h . IPettini et all ()2008[ ) suggested the 



large discrepancy in high-redshift DLAs could be due 
to delaye d enrichment of N with respect to O, esti- 
mated by iHenrv et ahl ()2000D to be --250 Myr, which 
is a more considerable fraction of time at high-redshift 
(see IPettini et al.) ()2008D for a more detailed discussion). 
Interestingly, all of our low-redshift systems also lie to 
left the transition point of secondary production seen in 



the H II regions of local spiral galaxies. The large spread 
in [N/a] is similar to that seen in high-redshift systems. 
This result seems puzzling if we expect these systems to 
have had sufficient time to begin secondary production. 
Despite this, the a abundances of the low-redshift DLAs 
and sub-DLAs are higher than that of the high-redshift 
systems, as one would expect. 

6. AFFILIATED GALAXIES 

The lack of information about the host galaxies of 
DLAs and sub-DLAs has been an obstacle that has lim- 
ited the exploitation of these systems for probing galaxy 
evolution. The deficiency in the number of low-redshift 
systems discovered to date has only hampered these ef- 
forts more. Some progress has been made recently us- 
ing photometric redshifts and colors to identify affiliated 
galaxies for 80 Mg II sel ected absorption systems at red- 
shift 0.1 < Zabs < 1-0 (|Rao et al.ll2011|) . However, our 
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Figure 9. Velocity plots for several lines in the HIRES spectra of 
J1619+3342. The dashed lines indicate the velocity components, 
the data is presented in black, the error array is given in blue, and 
the best fit is brown. Rest velocity corresponds to z = 0.0964. 

sample offers a unique opportunity to identify affiliated 
galaxies for systems not biased on Mg II absorption. 

To date, deep imaging has only been obtained for one 
of our absorbers. Imaging and speetroscopy of the field 
around J1009+0713 using the HST Wide Field Cam- 
era 3 and Keck/LR IS, respectively, was presented in 
iMeiring et al.l ()2011l ). This study revealed that none of 
the nearby galaxies are at the redshift of the DLA. We are 
limited to using the SDSS images for the other 7 sources 
to identify affiliated galaxies. Fortunately, spectra of sev- 
eral of the neighboring galaxies have been obtain e d wit h 
Kcck/LRIS as part of a survey by iWerk et all (j2011[ ). 
giving us accurate redshifts for some of the galaxies in 
the fields of our QSO. The redshifts determined from this 
survey are accurate to within 30 km s~^. For six of the 
galaxies observed, met allicities we r e also estimated using 
calibrati on method s of 'McGaugh| ()199lD and/or the N2 
index of lPettini fc Pagel (,2004) . The properties of these 



galaxies are summarized in Table 1101 The emission spec- 
tra for these gala xies are reproduced in Figure [13] from 
IWerk et all (l20Tll) . 

Affiliated galaxies have been found for three of the sub- 
DLA systems studied in this work. SDSS images of each 
of our absorbers, excluding J1009+0713, are shown in 
Figure [HI with candidate host galaxies indicated with a 
red box. Galaxies are identified using the position angle 
in degrees from the QSO, indicated by the first number, 
and the projected separation in arcseconds (impact pa- 
rameter) from the QSO, indicated by the second number. 
The large impact parameters for all of these galaxies sug- 
gest that the sub-DLAs are likely to be in their outer ha- 
los. There are a large number of dimmer features which 
have yet to be observed spectroscopically that may be 
affiliated galaxies. 

Two galaxies in the field of J0925-(-4004, identified as 
193_25 and 196_22, were found to be located at a sim- 
ilar redshift to the sub-DLA, making them good candi- 
dates as host galaxies. Galaxies 193_25 and 196_22 have 
impact parameters of p = 92 and 81 kpc, respectively, 
with respect to the QSO. Galaxy 196_22 appears to be 
more massive than its neighbor and also more luminous, 
with a L > as determined from its (jBlanton et al.l 
I2003f) . The metaUicity was estir nated for galaxy 193_25 
using the calibration methods of iMcGaughl (|1991[ ) to be 
about [0/H]= 0.12 ± 0.15, which is ~ 0.4 dex higher 
than the sub-DLA metaUicity. It is well known that the 
metallicity in galaxies can vary from their morphological 
type and also from the dista nce from galactic center. A 
study by iChen et all (|2005[ ) found a radial gradient of 
—0.041 ± 0.012 dex kpc~^ from the galactic center to a 
30 kpc radius for an ense mble of six g a laxy- DLA 
pairs. However, recent work bv iPeroux et al.l (|2011[ ) has 
shown that there is a large scatter in radial gradients and 
also cases of positive gradients in some systems. The ra- 
dial gradient necessary to explain the observed metallic- 
ity decrement from 193_25 would be —0.004 dex kpc~^, 
an order of magnitude lower than the rate found for the 
DLAsinlChen et al. (2005). However, since the sub-DLA 
is at a considerably large impact parameter, applying a 
gradient seems inappropriate. Given the proximity of 
these galaxies, it is more likely this gas was stripped from 
one of these systems during an interaction. 

Three galaxies in the field of J0928-I-6025, identified 
as 110_35, 129_19 and 187_15, were found to be located 
at a similar redshift to the sub-DLA, making them good 
candidates as host galaxies. Galaxies 110_35, 129_19 and 
187_15 have impact parameters of p = 89, 48 and 38 
kpc, respectively, with respect to the QSO. It is worth 
noting that these three galaxies are part of a larger clus- 
ter centered about 750 kpc away. Galaxy 110_35 ap- 
pears to be the most massive and luminous of the three. 
The metallicity of galaxy 129 19 was measured using 
th e method of [N IcGaugB (j 19911 ) and also the N2 index 
of IPettini fc Pagel (|2nol to be [0/H]= -0.22 ± 0.15 
and [0/H]= 0.07 ± 0.35, respectively. These values are 
below the metallicity estimated for the sub-DLA us- 
ing the corrected abundance of Fe. The metallicity of 
galaxy 187_15, the closest t o the sub-DLA, was m ea- 
sured using the N2 index of IPettini &: Pagell (|2004[ ) to 
be [0/H]= —0.38 ± 0.35, which is also lower than the 
sub-DLA metallicity. The complicated kinematic strue- 
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Figure 10. Velocity plots for several lines in the COS spectra of J1619+3342. The dashed lines indicate the velocity components, the data 
is presented in black, the error array is given in blue, and the best fit is browrn. Unrelated features or blends are indicated with magenta 
circles. Fe II A 1143.23 is blended with SUA 1253.81 in the MW. Rest velocity corresponds to ^ = 0.0964. 
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ture of this sub-DLA (required 13 veloeity components 
to fit) may be the result of an interaction between these 
galaxies. If this is the case, the sub-DLA could be tidal 
debris from one of these galaxies. Observational support 
for such a scenario has be e n seen in deep HST obser- 
vations by iKacprzak et al.l (j2010[ ). who observed large 
stellar streams and disturbed morphologies in a group 
of galaxies associated with a DLA, suggesting that the 
absorption system arises from tidal debris in the group 
environment. 

A single galaxy in the field of J1435-I-3604, identified 
as 68_12, was found to be located at a similar redshift 
to the sub-DLA, making it a good candidate as a host 
galaxy. Galaxy 68_12 has an impact parameters of p =38 
kpc with respect to the QSO. Interestingly, this galaxy 
shows a high star formation rate (SFR'^19 M0 yr~^) and 
also appears to be very luminous, with a L > L^, as de- 
termined from its M^- The metallicit y of galaxy 68 _1 2 
was measured using the methods of iMcGaughi (flofl 
and iPettini fc Pagell (I200I to be [0/H]= 0.12 ± 0.15 
and [0/II]= 0.08 ± 0.35, respectively. These values are 
about 0.5 dex higher than the metallicity we found for 
our sub-DLA. The radial gradient necessary to explain 
this observed metallicity decrement would be —0.013 dex 
kpc~^, or about o ne-third the rate found for the DLAs in 
iChen et al.l (|2005D . If 68_12 is the host, the distance from 
galactic center may explain the metallicity decrement. 
However, applying a gradient may not be appropriate at 
such large impact parameters. 

It is quite interesting that we find the potential hosts to 
these absorbers at such large impact parameters with the 
QSO. However, the fact that our candidates so far are at 
a relatively large impact parameter is not surprising for 
three reasons: (1) geometrically, it seems more likely to 
find an absorber intersecting the halo of a galaxy because 
of the relative size of the halo with respect to the disk 
or central region; (2) our sample was selected to be UV 
bright and since the QSO light will dim from extinction 
by the larger amounts of dust in a disk of a galaxy, it is 
possible those passing through a disk could be missed; (3) 
galaxies which arc at very small impact parameters from 
the QSO will not show up in the SDSS images because 
of the brightness of the QSO (however, small-impact pa- 
rameter galaxies can be identified by "composite" SDSS 
spectra showing emission from both the background QS O 
and the foreground galaxy, e.g.. lBorthakur et al.l (|2010( )). 
Deeper observations will need to be obtained to deter- 
mine if there are galaxies close to the QSO sightline. It is 
also worth noting that log TVhi was fou nd to be inversely 
correlated with p at the 3ct level by iRao et ahl (|2011[ ) 
and this seems consistent with only our sub-DLAs being 
found at large impact parameters. 

From the SDSS images it seems reasonable to conclude 
that none of the DLAs are affiliated with L, galaxies. An 
galaxy with a diameter of 30 kpc placed at a redshift 
of Zabs = 0.1 would be very bright, = 17.1, and have 
sufficient angular separation (13.4") to be distinguishable 
from the large PSF of the QSO in the SDSS images. It 
is likely that these systems are associated with dimmer 
galaxies in our field that have yet to be spectroscopi- 
cally observed. It is also important to perform deeper 
imaging of all of our fields, similar to that obtained for 
J1009-I-0713, to identify possible low-luminosity candi- 



dates at small impact parameters to the QSO. 

7. SUMMARY 

As part of our survey of galaxy halos with the Cos- 
mic Origins Spectrograph, 5 sub-DLAs and 3 DLA ab- 
sorption systems were serendipitously discovered. This 
sample nearly doubles the current sample of low-redshift 
QSO absorption line systems discovered to date. Accu- 
rate column densities and abundances were determined 
using a profile fitting procedure because of the low S/N 
of some of our sources and also to reduce the effects 
of the broad LSF of the COS spectrograph. The kine- 
matic structure of most of the absorbers appeared to 
be relatively simple, requiring only a few of components 
over a short velocity spread. Only J0928-f0625 and 
JlOOl-f 5944 showed complex kinematics, requiring more 
components spread over a larger velocity range. How- 
ever, it is likely that only the Keck/HIRES spectra re- 
flect the true component structure of the absorber and 
so further observations are needed to perform kinematic 
analysis of all of these systems. 

The metallicities of our observed systems show a large 
scatter in metallicity from near solar to one-tenth so- 
lar. Ionization corrections were made for the sub-DLAs 
using the CL OUDY ionization code. However, it has 
been shown by iMilutinovic et all ()2010f ) that the simpli- 
fication of a single phase medium can lead to a system- 
atic error in the result, in which the ionization correction 
factors could be too small. This in turn would lead to 
corrected abundances that may be too high. When com- 
bined with other sources from the literature, we see a 
trend of low-redshift DLAs being underabundant to so- 
lar metallicities. Sub-DLAs systems tend to have higher 
metallicities on average than DLAs at all redshifts. Un- 
fortunately, the sampling of absorptions systems are still 
highly skewed toward higher redshifts. With the current 
sample size, only 10% of DLAs and 23% of sub-DLAs are 
at look-back times less than half the age of the Universe. 
For look-back times less than one-quarter of the age of 
the Universe, this drops to just 3% of DLAs and 11% 
of sub-DLAs. A more intuitive approach to analyze the 
chemical evolution of the Universe would be to divide up 
the observed DLAs and sub-DLAs into equal-sized age 
bins. The number of low-redshift absorptions systems 
must increase if we want to make any substantial claims 
as to the nature of chemical evolution. Fortunately, this 
task is slowly taking shape with the new COS instrument 
aboard HST. 

For each of our low-redshift absorbers, except perhaps 
J1616-I-4154, the amount of N observed was significantly 
underabundant relative to a-process elements. Analy- 
sis of the [N/a] ratio for these systems showed that the 
majority lie below the transition point of secondary pro- 
duction levels seen in H II regions of local spiral galaxies. 
If these absorbers are associated with the disks of spiral 
galaxies then this result seems puzzling. It is expected 
that such low-redshift systems would have adequate time 
for secondary production levels to begin. 

Observations of low redshift DLAs and sub-DLAs are 
proving fruitful in understanding the nature of the host 
galaxies. Of the 8 systems we observed, we were able to 
identify affiliated galaxies for 3 sub- DLAs in the SDS S 
images using spectra from a survey bv lWerk et al.l ()201lD . 
The large impact parameters of these galaxies to the 
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Figure 11. [Z/H] vs look-back time for our low-redshiift systems along witii those from tiie literature. The red circles and green triangles 
correspond to unweighted mean metallicity of each bin for the DLAs and sub-DLAs, respectively. Horizontal bars denote the range in 
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Table 10 

Properties of possible host galaxies of sub-DLAs 



Field 


ID" 


, h 




At)'' 


SFR 


[O/H]- 


Log(M.) 


Mr 








(kpc) 


(km s-i) 


B aimer 


[Oil] 


M91-^ 


PP049 






J0925-I-4004 


193.25 


0.2467 


92 


-300 


0.86±0.15 


0.56±0.16 


0.12±0.15 




10.39 


-20.434±0.021 




196.22 


0.2475 


81 


-60 


<0.57 


<0.45 






11.07 


-21.546±0.013 


J0928+6025 


110.35 


0.1540 


89 


60 


<0.03 


<0.04 






10.56 


-20.411±0.012 




129.19 


0.1542 


48 


120 


2.89±0.37 


4.67±1.35 


-0.22±0.15 


0.07±0.35 


9.86 


-19.728±0.021 




187.15 


0.1537 


38 


-30 


0.45±0.05 


0.31±0.09 




-0.38±0.35 


9.33 


-18.664±0.039 


J1435+3604 


68.12 


0.2024 


38 


-60 


18.96±2.28 


15.57±4.43 


0.12±0.15 


0.08±0.35 


10.87 


-21.461±0.011 



Note. Measurements taken from lWerk et all 1)20111 ) 

" Galaxy Identifier, where the first number is the position angle in degrees from the QSO and the second number is the projected separation 
in arcseconds (impact parameter) from the QSO. 
Uncertainty in redshift values is 30 km s~^. 

Projected separation between the galaxy and QSO in kpc, calculated in the rest frame of the galaxy. 

Relative velocity difference of galaxy with respect to Zabs- Uncertainty is the same as redshift (±30 km s^^). 
<= For reference [O/H]0 = 3.31. 
■' Abundance using [McGaughl ||1991|) calibration. 
^ Abundance using the N2 index of lPettini fc Pagell l[2003 '). 



sub-DLAs suggest the absorbers lie in their outskirts. 
None of the DLAs appear associated with galaxies, as 
these would be easily observed in SDSS images. A large 
number of dimmer objects in our fields have yet to be 
measured spectroscopically and could likely yield more 
galaxy affiliations. More observations will be needed 
to properly identify the host galaxy of each of our ab- 
sorption systems. Follow-up studies will ultimately im- 
prove our ability to understand the absorber-galaxy rela- 
tionship and advance theories on chemical evolution and 



galaxy formation. 
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Figure 13. The ID reduced, flux-calibrated spectra for the galaxies included in this work. We represent the dichroic with a shaded area 
near the observed wavelength 5000 A. Observations from lWerk et all 1)20111 ). 
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APPENDIX 
COLUMN DENSITY MEASUREMENTS 

Here we breifly describe our method for using profile fits on each of the absorption systems. The velocity components 
and b values for each system, as well as the column densities of detected species at each of these components, are shown 
in Table HH 

The Sub-DLA m SDSS J0925+4004 

The velocity centroids and b values for this system were determined using O I because we detect five lines at differing 
line strengths for this species. The detection of so many O I lines is often difficult because the majority of its lines 
lie below 1000 A, which is typically outside instrument coverage unless the system is significantly redshifted. The 
kinematics of this system appear simple with only 2 components needed to properly fit the observed profile at COS 
resolution. The best fit parameters from O I were then used to fit N I, Si II and Fe II, and the results can be seen 
in Figure [2] This system is unique among our sample as it shows numerous molecular hydrogen absorption features. 
These absorption lines will be analyzed in a future paper. 

The Sub-DLA m SDSS J 0928+ 6025 

Optical spectra of J0928+6025 were obtained with the Keck/HIRES spectrograph. The absorption lines detected 
by Keck for this system arc shown in Figure [3l The velocity centroids and b values for this system were determined 
using Mg I. The kinematic structure of this system is quite complex and required 9 components to properly fit the 
observed profile. An additional 4 components were necessary to fit weaker outer components appearing in the much 
stronger Mg II lines. The Ca II line was independently fit, as it shows slightly different kinematics than Mg I, and 
the result is presented in Table [121 Since the resolution of the COS instrument is not high enough to distinguish all 
the components found in Mg I, the attempted COS fits utilized only the components at v = —39.52,-18.95, 3.5, 36.51 
km s~^, as these appear to be the strongest features in the Keck data. Unfortunately, the lines observed in the COS 
data were too strongly saturated or were too noisy for reasonable fits to be made of the observed profile. This can be 
seen in Figure Therefore, apparent column densities are adopted for each line in the COS data. For this case, lines 
were determined to be unsaturated if the apparent column densities of two or more lines with different line strengths 
were consistent. 

The Sub-DLA in SDSS J1001-h59U 

The velocity centroids and b values for this system were determined using O I because of the large number of observed 
lines at differing line strengths. The kinematics of this system appear complex with 4 components over a large velocity 
spread being needed to properly fit the observed profile at COS resolution. The O I parameters were then used to fit 
N I, Si II, PII, and Fe II, and the results can be seen in Figure [SJ This sub-DLA shows a Lyman limit break in the 
continuum at a wavelength of 1195 A in the COS rest frame. 

The Sub-DLA m SDSS J1435-h3604 

The velocity centroids and b values for this system were determined using N I because there are 6 observed lines 
and most of the O I lines appear saturated or blended. The kinematics of this system appear simple with only 2 
components needed to properly fit the observed profile at COS resolution. The velocity parameters from the N I fit 
were then used to fit Si II, S II, and Fe II, and the results can be seen in Figure [6l 

The Sub-DLA m SDSS J1553+3548 

The velocity centroids and b values for this system were determined using Si II since it was the most extensively 
covered line. The kinematics of this system appear very simple with only a single component needed to properly fit 
the observed profile at COS resolution. The Si II parameters were then used to fit N II and Fc II, and the results can 
be seen in Figure [T] 

The DLA in SDSS J 161 6-1-4 154 

The velocity centroids and b values for this system were determined using Fe II. The kinematics of this system appear 
simple with only 3 components needed to properly fit the observed profile at COS resolution. The parameters from 
the Fe II fit were then used to fit C II*, P II, S II, and Fe II and the results can be seen in Figure [51 

The DLA in SDSS J 1619+3342 

Optical spectra of J1619+3342 were obtained with the Keck/HIRES spectrograph. The absorption lines detected 
by Keck for this system are shown in Figure |9l The Mg II lines lie at the edge of the detector and so have very low 
S/N. The velocity centroids and b values for this system were determined using Mg I. The kinematic structure of this 
system is simple and required only 2 components to properly fit the observed profile. For very strong transitions, such 
as Mg II, C II, and Si HI, an additional component at w = —143.0 km s"^ seems to appear. However, this additional 
component was not used in the fits of the weaker lines. The Mg I fit parameters were then used to fit N I, Si II, P II, 
S II, Ca II, Fe II, and Fe HI and the results can be seen in Figure [TUj Since the resolution of the COS instrument is 
lower than that of Keck, some of the COS fits utilize only the component &i v — —21.0 km s^^, since it appears to 
contain most of the column density. 



Table 11 

Component structure and column densities for each absorber 



QSO 


Vel 
(km s-i) 


b 

(km s-i) 


logWcii* 
(cm-2) 


logWNI 

(cm-2) 


logTVNii 
(cm-2) 


logA^oi 
(cm-2) 


logTVMgi 
(cm-2) 


logAfsiii 
(cm-2) 


logWpii 
(cm-2) 


logNsu 
(cm- 2) 


logWcall 

(cm-2) 


logAfxiii 
(cm- 2) 


logTVpell 

(cm-2) 


logWFelll 

(cm-2) 


J0925+4004 


-34±4 
41±2 


24±4 
20±3 




14.14±0.06 
14.63±0.05 




15.36±0.06 
15.82±0.11 




14.32±0.07 
14.32±0.09 










13.16±0.76 
14.04±0.11 




J0928+6025 


-39.5±0.6 
-29.6±0.5 
-24.4±1.1 
-19.0±2.0 
-13.5±3.6 
3.5±1.0 
14.5±1.0 
28.1±9.4 
36.5±0.5 


6.3±0.9 
0.5±0.5 
2.1±1.7 
2.0±1.9 
4.0±3.1 
8.7±2.1 
2.1±1.6 
8.5±7.2 
3.8±1.1 










12.00±0.05 
10.96±2.42 
11.53±0.23 
11.32±0.98 
11.56±0.49 
12.15±0.07 
11.25±0.37 
11.59±0.66 
11.98±0.23 
















JlOOl+5944 


-38±1 
20±2 

63±18 
92±6 


20±2 
18±4 
13±13 
22±9 




13.44±0.09 
13.06±0.19 

12.50±0.71 




15.44±0.02 
15.15±0.03 
13.86±0.14 
14.01±0.11 




14.50±0.04 
14.26±0.04 
13.08±0.07 
13.43±0.04 


12.75±0.07 
11.93±0.33 








14.05±0.05 
13.88±0.06 
13.07±0.33 


13.88±0.08 
13.62±0.14 
13.17±0.34 


J1435+3604 


-31±2 
-4±1 


9±4 
5±1 




13.58±0.08 
14.56±0.16 












14.00±0.23 
14.47±0.15 






13.54±0.20 
14.36±0.12 




J1553+3548 


-29±1 


29±4 






14.16±0.07 






14.22±0.05 










14.01±0.07 




J1616+4154 


-106±20 
-59±10 
10±2 


20±15 
28±16 
24±2 


13.43±0.17 
13.79±0.09 












12.39±0.98 
13.42±0.10 


14.57±0.47 
15.29±0.12 






13.79±0.61 
14.32±0.23 
14.89±0.03 


14.24±0.06 
14.18±0.06 


J1619+3342 


-29.8±2.5 
-21.0±1.4 


4.4±2.1 
4.1±1.2 




14.63±0.12 






11.89±0.26 
12.24±0.18 




13.17±0.19 


15.08±0.09 


11.55±0.04 
12.36±0.01 


11.11±0.19 
11.83±0.04 


14.25±0.16 


13.45±0.37 
13.78±0.28 



Note. This table is just component-by-component column densities from the profile fit and does not include AOD column density estimates or limits. 
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Table 12 

Parameter fit to Ca II in the =0.1538 SDSS J0928+6025 sub-DLA 



Vel 


b 


logAfcaii 


(km s^-'-) 


(km s-i) 


(cm~^) 


-39.8±0.3 


4.8±0.5 


11.98±0.03 


-29.4±0.9 


3.2±1.6 


11.38±0.20 


-21.5±0.5 


3.5±1.0 


11.80±0.11 


-12.5±0.7 


4.5±1.3 


11.75±0.11 


0.2±2.2 


4.7±1.9 


11.90±0.32 


8.4±2.6 


5.2±2.6 


11.92±0.32 


18.2±0.8 


2.4±1.3 


11.26±0.24 


31.7±2.1 


8.9±1.7 


11.89±0.12 


36.4±0.2 


3.3±0.4 


12.18±0.06 



